Self-consistent calculation of the autolocalization barrier for quasiparticles in 

anisotropic crystal. 



A.Shelkan 

Institute of Physics, University of Tartu, Riia 1^2, EE2400 Tartu, Estonia 

(February 7, 2008) 

The energy of the electron wave packet interacting with lattice distortion, is considered in 
anisotropic crystal. Anisotropy of the electron and phonon spectra as well as of the electron-phonon 
interaction are taken into account. The height of the barrier between free and self-trapped states 
is calculated in dependence on the anisotropy parameters. The calculation is done numerically, 
using continual aproximation. The analytical solution is obtained for some cases of quasi-two and 
quasi-one-dimensional spectra. 
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There is a general agreement that in systems with distortion-type electron-phonon coupling, an electron (hole, 
exciton) forms the self-trapped l^nolaronic-type) state. Usually, this state is separated from the fr.ee one by an 
autolocalization barrier (AB , seeon). Tlie theory of the AB has been developed far isotropic systemscl. During the 
last years polarons in anisotropic systemsa, especially in superconducting materialsa, are widely discussed. Therefore 
it is of interest to consider how the lattice anisotropy influences the AB . Briefly this problem was considered ino, 
where the one- and two- dimensional limits were studied. The phonon spectrum and electron-phonon coupling have 
been taken to be isotropic. It was shown, that within these limits the height of the barrier E ~ {jnim^mz) , i.e 
barrier vanishes in one- and two-dimensional cases. However, in three-dimensional crystals with finite anisotropy the 
barrier has not been studied. n 

In this paper, we apply the adiabatic theory of AB, developed irfl, to anisotropic crystals. To make the problem 
easier we consider the crystal with symmetry not lower than romboedrical (i.e. with diagonal tensors of effective mass 
and elastic constant). To simplify the operator of the electron-phonon interaction we suppose the localization to be 
connected with the displacement of the particles along only one direction (signed z), that corresponds to the direction 
of the lowest dispersion of electron zone. This approximation is justified by the fact that in the wide zone limit 
considered here, the influence of the electron-phonon interaction on the energy of the electron is inverse proportional 
to the width of the zone. Besides, we suppose that the energy of the AB essentially exceeds the average energy of 
acoustical phonons and the radius of the barrier state is much larger than the lattice constant. These assumptions 
allow us to use the adiabatic approximation and to solve the problem in the continual limit. 

The energy functional readscl: £ — £ e + £ p + V 
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Here C(r) denotes the wave function, and £ e - the cinetic energy of the electron, £ p is the lattice elastic energy, V - 
the energy of the electron-phonon coupling , and u(r) labels the component of the atomic displacement. In this model 
exists a self-consistent state that corresponds to the barrier and that is the saddle point of the energy functional . To 
determine this point we vary the energy functional initially with respect to u(r). Then we choose the wave function 
C(r) in exponential form 

W^'ClW'^-^l, (2) 

and consider R(i) as variational parameters, corresponding to the radius of the electron state. As a result we get for 
the energy functional 



1 



3 \2 



=1 R i 



where 
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For reducing the number of free parameters we simplify the problem: we restrict ourselves to a system that is 
isotropic in the (x, y) plane, i.e. we suppose that R\ = R% = R; J\ = J2 = J] @i= fh = (3- Parameters (3 correspond 
to the elastic constants^: (3 = X ZX zx ; 0a = Xzzzz ■ In the isotropic systems \ ZZ zz - X zzxx - 2X ZXZX = 0. In the 
approximation of central forces \ ZX zx = X ZZXXl and the elastic isotropy condition takes the form /3 3 = 3/3. Evidently 
the condition of spatial isotropy for the electron spectrum is J 3 = 3 J . It is convenient to characterize the violation 
of spatial isotropy by the dimensionless parameters 7 = J/J3 , a — 3/3/ fy. 

Then 

£ = J 3 "52 + ^2 " TT7T / 7~~r? \ 7 • ( 5 ) 
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Let us turn to the new variational parameters d and 0, introduced according to: R 2 = dsm.8,R^ = dcosO. After 
variation with repect to d the energy functional takes the form: 

8 = ^X^ctge - ^ + ■ (6) 
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Now the problem is reduced to the minimization of the energy with respect to one single parameter 9 - this minimal 
value E is the height of the barrier. We perform this minimization numerically. 

From (6) one can see that, besides a and 7 , the result depends also on J, (3 and A. Our goal is to elucidate 
how the anisotropy parameters influence the height of the barrier. For this purpose we shall look for the ratio 
E/Eq of the heights in anisotropic and isotropic cases E/Eq. For the scale parameters it is convenient to choose 
,P S = J™ = (2 J + J 3 )/3, f3 ls = I3\ s = (2/3 + #0/3, X ls = A. Then the ratio E/E Q depends only on the anisotropy 
parameters a and 7 . Similarly we introduce the ratios R/Rq and R3/R0. The results of numerical calculations of 
the ratios E/Eq as functions of the anisotropy parameters are presented in Figs 1-2. One can see that the maximal 
height of AB corresponds to 7 ~ 1 and a ~ 0.1. In the case of the quasi-one-dimensional (7 << 1) electron spectrum 
barrier vanishes slower than in the quasi- two-dimensional limit (7 >> 1) . The transition to the one-dimensional 
phonon spectrum (a << 1) only weakly affects the height of the barrier, and for large 7 doesn't affect the barrier at 
all. 

We note that for the parameters of anisotropy being used in our numerical calculations, the ratios R/Rq and R-j/Rq 
are of the order of one or larger. Therefore the continual approximation yields a correct description of the AB in 
anisotropic crystals if it is applicable in corresponding isotropic crystals. For a >> 7 and 7^0 the ratio R/Rq —> 
, i.e. in this limit we cannot use the continual approximation. Besides, in anisotropic cases, i.e. when AB vanishes, 
then either one or both radii go to infinity. 

In the limiting cases of large and small values of the ratio a/7|£>ne can obtain analytic solutions. Denoting aR 2 /3R 2 
as u, the integral in (5) as f(u) and using formula 3.197.1 fromQ one obtains 

f{n) - / * ( . + ,4 + - 1 ~ l F 1 -") ■ < 7 » 

where F is the hypergeometric function. _ 
Let us consider first the case a << 7 and assume that u « 1. Transforming (7) according to 9.132.1 fromQ, 
neglecting the second order terms in it 1//2 and variating with respect to R and R3, we obtain 
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From (8) one can see, that u ~ a/7 . The condition u << 1 is fulfilled, so the corresponding assumption is justified. 

Now let us consider the case a » 7. Here we assume that u » 1. Using the properties of hypergeometric 
function we rewrite (7) as 

/(tO^F^Mj^l-u) (9) 

Neglecting the second order terms in 1/u and variating with respect to R and R3 we get 

E a 2 A a x ~ 2 



So (6a+ 1) 2 (2 7 + l) 3 V ln 7 
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One can see from (10) , that the condition u >> 1 is fulfilled. In the case of the quasi-two-dimensional phonon 
spectrum (a » 1), E/E goes to zero as (In a) -2 , and R/Rq and R3/R0 goes to infinity correspondingly as (In a) 3 / 2 
and In a . The transition to the quasi-one-dimcnsional phonon spectrum (a << 1) slightly influences the height of 
the AB. If the electron spectrum is quasi-two-dimensional (7 >> 1), then E/Eq goes to zero as 1/7, and R/Ro goes 
to infinity as 7 1 / 2 . This means, that the effect of anisotropy of the electronic spectrum to the AB is much stronger, 
than of the anisotropy of the phonon spectrum. 

Note that within the framework of this model, the consideration of the limit, which correponds to the one- 
dimensional electron spectrum (7 << 1) is not justified. Indeed, when 7 goes to zero, then as one can see from 
(12), the ratio R/Ro also goes to zero and the continual approximation becomes inapplicable. 

In conclusion, in this communication the influence of the anisotropy of electron and phonon spectra and also of 
electron-phonon interaction on the height of the barrier of a phonon polaron is investigated. It is shown, that the 
anisotropy of the electron spectrum stronger affects the AB, than the anisotropy of the phonon spectrum. In the 
quasi-two-dimensional limit the height of the barrier goes to zero. We emphasize that the disapearence of the barrier 
in the quasi-two-dimensional case is characteristic for the phonon polaron. On the contrary, as it was pointed out inta, 
a spin-polaron model characterized is by the a finite barrier also in a two-dimensional AF lattice. The reason of this 
difference between phonon and spin polarons lies in the strongly nonlinear nature of the localized spin excitation: one 
turned spin forms a stable topological defect in the AF lattice even without a charge carrier. Therefore a large-size 
wave-packet of a free particle with one turned spin represents a metastable state ( unlike the wave packet with a 
local lattice distortion, which is fully unstable already in the time scale of the order of the vibrational period ) . This 
difference in the properties of the phonon and spin polarons must be taken into account when discussing the transport 
properties of the high-T c superconductors. 

I wish to express my gratitude to D.Nevedrov for valuable discussions and stimulating of this publication. 
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Fig.l The ratio of the heights of the barrier in anisotropic and isotropic cases E/E vs anisotropy parameter of 
electron spectrum 7 for various anisotropy parameters of phonon spectrum a. 
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Fig. 2 The ratio of the heights of the barrier in anisotropic and isotropic cases E/E vs anisotropy parameter of 
phonon spectrum a for various anisotropy parameters of electron spectrum 7. 
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